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Abstract 


Equations of motion are derived for use in simulating a spacecraft or other com- 
plex electromechanical system amenable to idealization as a set of hinge-connected 
rigid bodies of tree topology, with rigid axisymmetric rotors and nonrigid append- 
ages attached to each rigid body in the set. In conjunction with a previously pub- 
lished report on finite-element appendage vibration equations, this report provides 
a complete minimum-dimension formulation suitable for generic programming for 
digital computer numerical integration. 


vi 


JPL TECHNICAL REPORT 32-1576 



Dynamic Analysis of a System of Hinge-Connected 
Rigid Bodies With Nonrigid Appendages 


I. Introduction M 

In a previously published report (Ref. 1), there appear equations of motion which 
characterize the small, time-varying deformations of an elastic appendage attached 
to a rigid body experiencing arbitrary motions in inertial space. The flexible append- 
age is modeled as a set of deformable elastic elements possessing distributed mass 
and interconnected at n nodes, with a rigid nodal hody appearing also at each node. 
This finite-element model has 6n degrees of freedom in deformation, corresponding 
to the degrees of freedom invested in the nodal bodies; deformations of the inter- 
nodal elastic elements are established by assigned interpolation functions. The 
purpose of Ref. 1 is to establish the structure of the 6n deformation equations, in 
order to permit consideration of coordinate transformations which might introduce 
the possibility of coordinate truncation and the consequent representation of elastic 
appendage deformations in terms of distributed or modal coordinates numbering 
much less than 6n. With this objective accomplished, Ref. 1 terminates, leaving a 
set of equations of motion which are incomplete in the sense that they are insuffi- 
cient to determine the kinematic variables characterizing the motion in inertial 
space of the rigid base to which the flexible appendage is attached. 

It is the purpose of the present report not only to complete the dynamic analysis 
begun in the earlier work, but to do so in a way that encompasses a wide class of 
vehicles, namely, those amenable to idealization as a set of n + 1 rigid bodies inter- 
connected by n line hinges (implying tree topology), with the possibility of rigid 
axisymmetric rotors and arbitrary nonrigid appendages attached to each rigid body 
in the set. 
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The results of Ref. 1 provide the vibratory deformation equations for each elastic 
appendage in the system, and the transformation to distributed coordinates which 
is appropriate in each case. The new scalar equations to be derived in this paper 
number at least n + 6, being descriptive of the inertial translations and rotations of 
one reference body and the relative rotations about the n line hinges, with an addi- 
tional equation being added for each axisymmetric rotor in the system. In this 
respect there is a strong parallel between the results of the present report and those 
obtained by extending the Hooker-Margulies equations (Ref. 2) as suggested by 
Hooker (Ref. 3) in order to eliminate unwanted kinematic constraint 1 torques; the 
significant difference is that Hooker and Margulies considered only point-connected 
rigid bodies in a topological tree, whereas in the present report the basic elements 
in the tree are substructures which are rigid in part but include rotors and arbitrary 
nonrigid appendages. 

A somewhat restricted version of these results is presented in Section VI of this 
report in a matrix form which has an affinity with the multiple-rigid-body formalism 
developed by Roberson and Wittenburg (Ref. 4) in chronological parallel with the 
derivation of the Hooker-Margulies equations (Ref. 2), although the form of the 
present equations lacks the aesthetic qualities of those in Ref. 4. 

A digital computer program for the numerical integration of the general equa- 
tions reported here is under development at the Jet Propulsion Laboratory by 
G. E. Fleischer. The result is to be a generic program, suitable for the dynamic 
simulation of a wide class of spacecraft. Many features of the following derivation 
have been adopted so as to minimize the labors of the users of this program. 


II. Mathematical Model 

Any problem of dynamic analysis must begin with the adoption of a mathe- 
matical model representing the physical system of interest. In what follows, it is 
assumed that the model consists of n + 1 rigid bodies (labeled / 0 , • ■ - , 4 n ) inter- 
connected by n line hinges (implying no closed loops and hence tree topology), 
with each body containing an arbitrary number (perhaps zero) of rigid rotors, each 
with an axis of symmetry fixed in the housing body, and moreover with the possi- 
bility of attaching to each of the n + 1 bodies a nonrigid appendage, with append- 
age attached to body 4%. The appendage itself can be modeled in a variety of 
ways without exceeding the scope of the final vector-dyadic equations in this 
paper; one might adopt a continuum model, a distributed-mass finite-element 
model, or a model admitting mass only in the form of nodal bodies or nodal 
particles. 

Specific choices of appendage model are made only in Section VI of this report, 
where an explicit set of matrix equations is presented in order to illustrate the 
transition from generic vector-dyadic equations to the more explicit matrix or 
scalar equations required for digital computer numerical integration. Although the 
equations derived here do not in their vector— dyadic form imply restricted append- 
age deformations, when specific appendage models are adopted it is assumed that 
their deformations relative to some reference state are “small” in the sense that 
terms above the first degree in deformation variables can be ignored. 
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L Kinematic constraint forces and torques are respectively those interaction forces and torques 
which maintain kinematic constraints, such as “no relative translation of two points” or “no devi- " 
ation from a prescribed relative rotation.” 
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If the actual connection between two massive portions of the physical system 
admits two (or three) degrees of freedom in rotation, then the analyst simply intro- 
duces one (or two) massless and dimensionless imaginary bodies into his model 
(as though they were massless gimbals). Since the number of equations to be 
derived here matches the number of degrees of freedom of the system, no price is 
paid in problem dimension by the introduction of imaginary bodies, and considera- 
ble simplification results in the user input format for the computer program. 

Each combination of a rigid body and its internal rotors and attached flexible 
appendage comprises a basic building block referred to here as a substructure ; 
thus, there are n 4- 1 substructures in the total system, so labeled that encom- 
passes and any rotors in 6 k . 

Figure 1 illustrates a typical spacecraft and its mathematical model; flexible 
appendages are cross-hatched. Note that the three physically distinct flexible struc- 
tures attached to the central body are all combined as the single flexible appendage 
of in the mathematical model, and the two-degree-of-freedom gimbal joint of 
the midcourse motor has been accommodated with the introduction of an imaginary 
massless rigid body Because the fifth substructure encompasses the flexible 
antenna as well as the rigid body <£ 5 , and because the mass properties to appear in 
the equations are those of the substructure (and not those of the rigid body), the 
mass properties of are irrelevant to the analysis. The masses of and the rotors 
within / 0 are also irrelevant, being absorbed in the mass of substructure <* 0 . 

Normally we must expect that the relative rotations between contiguous pairs of 
substructures will be governed by components of torque generated by electric 
motors subject to control laws. It is this expectation that underlies the transition 
from the physical system to the mathematical model in Fig. 1, This is not a uniquely 
satisfactory model, however, and if attention were to focus on vehicle attitude con- 
trol rather than on antenna pointing control, it might be desirable to replace the 
control systems governing the rotations between and and between and 
^0 by analytical expressions typical of elastic or viscoelastic connections. If this 
choice is made, it might be best to embrace the fourth and fifth substructures 
within the definition of the flexible appendage of <f 0 . This decision would result in 
a "better” set of modal coordinates (permitting simulation with fewer variables), 
as one can appreciate after careful study of Ref. 1. The importance of devising the 
most appropriate mathematical model for the task at hand deserves emphasis, since 
with the advent of automated derivation and integration of equations of motion 
this becomes the most important step in dynamic analysis. 


III. Derivation Procedure 

The derivation is described below. 

Step I. Isolate each substructure and apply 

F{, = QftyA' and T ? t = H ; ; = 0, 1, - * , ft (1) 

where, for the fth substructure, is the total resultant of all external forces, Q yij is 
the substructure mass, A* is the inertial acceleration of the substructure mass center 
Cj, T* t is the total moment resultant of all external forces referred to Cj , and H* is 
the inertial time derivative of the substructure angular momentum referred to Cy. 
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(a) VEHICLE 


ANTENNA- 



Fig. 1. A spacecraft and its mathematical model 
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Step 2. Combine the 6n +■ 6 scalar equations obtained in Step 1 so as to obtain 
n + 6 scalar equations which do not involve redundant variables or those substruc- 
ture interaction forces or torques which serve to maintain kinematic constraints. 

Step 3. Apply T = H to each rotor of the system (where symbol definitions follow 
naturally from Eq. 1), and dot-multiply each equation by a unit vector parallel to 
the symmetry axis of the corresponding rotor; the result is a set of scalar equations 
matching in number the rigid, axisymmetric rotors in the system. 

Step 4. Record the appendage deformation equations from Ref. 1 (or an alterna- 
tive source, depending on the appendage model), and, if a finite-element model is 
used, substitute the deformation coordinate transformations from Ref. 1 wherever 
the deformation variables appear in the preceding equations. Maximum generality 
is retained in this report by postponing such substitutions until the equations are 
in otherwise final form, suitable for any appendage model. 

Step 5. Specify all control laws in the form of ordinary differential equations in 
time, with control torque magnitudes or their equivalent as dependent variables. 


IV. Definitions and Notations 

Definitions and notational conventions are as follows (see Fig. 2): 

Def. 1 . Let n be the number of hinges interconnecting a set of n + 1 substructures. 



Fig. 2. Definitions for the Jcth substructure, with / < k 
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Def. 2 . Define the integer set £3 — {0, 1, * * 1 , n}. 

A 

Def. 3 . Define the integer set & = {1, * * * ? n}. 

Def. 4. Let <£ 0 be a label assigned to one rigid body chosen arbitrarily as a refer- 
ence body, and let £ u * * * , 4 n be labels assigned to the rest of the rigid 
bodies in such a way that if 4j is located between 4 0 and 4 k then 0 < / < k. 


Def . 5. Define dextral, orthogonal sets of unit vectors b*, b*, b* so as to be im- 
bedded in 4 k for k e £3, and such that in some arbitrarily selected nominal 
configuration of the total system, b* — b^for a ~ 1, 2, 3 and fc, /e£3. 


Def. 6. Define 



ke£B 


Def , 7. Define {i} as a column array of inertially fixed, dextral, orthogonal unit 
vectors i 1? i 2 , i 3 * 


Def. 8. Let C be the direction cosine matrix defined by 

{b°}-C{i} 


Def. 9 . Let go 0 = {b 0 } 7 ® 0 be the inertial angular velocity vector of 4 0y so that o>° 
is the corresponding 3X1 matrix in basis {b 0 }. 

Def. 10. Let c k be the mass center of the Jcth substructure, k e£3, 

Def . 11. Let /i k be a point on the hinge axis common to 4 k and 4 f for ; < k and 

k*g>. 

Def . 12. Let p kj be the position vector of the hinge point connecting 4$ and 4 k 
from the point o* occupied by c k when the kth substructure is in its nomi- 
nal state. 


Def. 13. Let c* be the position vector from c k to o k . 

Def. 14. Let p ft be the position vector to c k from the system mass center CM. 

Def. 15. Let X be the position vector to CM from an inertially fixed point ^ and 
let X = X> {i}. 

Def. 16. Let Q 7l k be the mass of the fcth substructure, for k e£3. 

Def , 17. Let (p)* be a generic position vector from o* to any point in the kth 
substructure. 

Def. 18. Let Q k be a point common to rigid body 4 k and flexible appendage a k . 

Def. 19. Let R* = {b be the position vector fixed in 4 h locating Q k with 
respect to o k . 

Def. 20. Let (r) fc = (b fc } r (f)* be a generic symbol such that R fe + (r) fc locates a 
typical field point in a k with respect to o k when the flexible appendage is 
in some nominal state (perhaps undeformed). For a discretized appendage 
let (r a )* = (b*} r ( r *) k locate the sth node in the nominal state. 
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Def . 21. Define the generic deformation vector (u)* in such a way that 2 

(p) fc ^ R* 4- (r) fc + (u) fc 

and 

( p) k = R k 4“ ( r) k 4- ( u) k 

For a discretized appendage a ky let (u*) fc = {b h } T (u*)* be the deformation 
vector for node s. 

Def. 22. Let g* = (b^g* be a unit vector parallel to the hinge axis through 

Def . 23. For ke£P, let y * be the angle of a g k rotation of with respect to the body 
attached at / % k . Let y* be zero when b* = b* (a = 1, 2, 3; ket 3). 

Def. 24. Let J* = {b k } T J k {b fc } be the inertia dyadic of the fcth substructure for o ky 
so that J* is time- variable by virtue of deformations. 

Def. 25. Let F fc = (b fc } r F k be the resultant vector of all forces applied to the fcth 
substructure except for those due to interbody forces transmitted at hinge 
connections. 

Def . 26. Let T k = (b*} r T k be the resultant moment vector with respect to c k of all 
forces applied to the fcth substructure except for those due to interbody 
forces transmitted at hinge connections. 

Def . 27. Let r* be the scalar magnitude of the torque component applied to in 
the direction of g k by the body attached at /i k . 

Def. 28. Let F = 2F k = { b°} T F be the external force resultant for the total system. 

ke£& 

Def. 29. Define the scalar such that for ke£B and s e & 

1 if /i s lies between <f 0 and 6 k 
0 otherwise 

(The n (n 4- 1) scalars e sk are called path elements.) 

Def. 30. Define Oft = 2 0 7lk> the total system mass. 

kctB 

Def. 31. Let C rj be the direction cosine matrix defined by {b r } = C ri {h i } 7 r i j€& B. 
(Note that in the nominal state, C rj = U , the unit matrix.) 

Def. 32. Let N kT denote the index of the body attached to <£& and on the path lead- 
ing to 4 r , and let N kk = k. (These are the network elements.) For notational 
simplicity, use N k for N kQ . 

Def. 33. For 3 re£Q ~ k, let L fcr = and let V th = 0. 

Def. 34. Define D» = -2 for k e £8. 

ft® 

Def. 35. Let b k be a point fixed in <4 such that D kk is the position vector of o k with 
respect to b k . (This point b^ is called the barycenter of the /cth substructure 
in the nominal state.) 



Superscripts on generic symbols such as p, r, and u will be omitted when obvious, as when the 
symbol appears within an integrand of a definite integral. 

3 For notational brevity, the set £8 — {&} is designated £8 — fc. 
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Def. 36. Define {b*} r D** = = D** + h k > for k,je(B. 

Def . 37, Define the dyadic 

K* = 2 Qftr (D fcr • D* r U - D*D*) 

where U is the unit dyadic, and define the corresponding matrix K* == 
{b*} • K* * { b k } T . 

De/. 38, Define 

•** = K* + J* and = {b*} • $** ■ {b*} T 

D#/. 39. Define 

<!>*> = -Q77 (d^-D^U - D /k D*') 

with 

{V} •**'• {b h y = -0% (C^D^C^Z)*' - D ik D k ’ T ) 

Def. 40. Let the system of forces applied to 4 k by the attached body 4} be equiva- 
lent to a resultant force F* passing through the labeled point {ft$ or fn) 
common to 4 k and 4 h plus a torque T*h 

Def . 41, Let t*> be the kinematical constraint torque applied to 4 k by 4 h in such 
a way that, with Defs. 40, 22, and 27, 

T*' = t*' + 8 iyfc T*g*-8*y,Tj g' 

where 

$jj = 1 and 8/ ; = Ofor t =4=/. 

Def. 42, Let ai* = {b h } T o> k be the inertial angular velocity of 4 k . 

Def. 43. Let £3 r be the rth neighbor set for r e £B, such that A: e £Q r if 4 k is attached 
to 4 r . 

Def . 44. Let <©,-* be the branch set of integers r such that r €fB jk if fe = N jr . Thus 
consists of the indices of those bodies attached to 4, on a branch 
which begins with 4 k . 

Def. 45. Let h* be the contribution of rotors in 4 k to the angular momentum of the 
fcth substructure relative to 4 k with respect to and let h* = h k * {b fc } . 


V. Derivations of Vector-Dyadic Equations 

In terms of the indicated definitions, Eq. (1) provides, for the rth substructure 

(re&)> 


F r +S f"-OMr(X + p') = 0 

9*&r 


(2) 


and, for the fcth substructure (fee £3), 

T* + 2 T*' + 2 (p* + c*) X W - H* = 0 (3) 

ft®* 

Here a dot over a vector implies time differentiation in an inertial frame of reference. 
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As shown in Ref. 5, pp. 32-33, the fcth substructure angular momentum H fe can be 
expanded from basic definitions as 



where superscripts k are dropped from c and p within integrals ranging over the 
substructure Here the open circle over a vector implies time differentiation in 
the reference frame established by 4- The contribution of a rotor to the integral 
last appearing in Eq. (4) may be recognized as a vector fixed in 4 k and parallel to 
the rotor axis, with a magnitude equal to the product of rotor spin axis inertia and 
the relative angular speed of the rotor relative to 4 k ; the vector sum of the contribu- 
tions to this integral of all the rotors in 4 k is here designated h*. Since p h = 0 for any 
field point in Eq. (4) can be written 


H fc “ J* • 4- X + h fc + / pXp dm (5) 

limiting the range of integration to the appendage Substitution of the inertial 
time derivative of Eq* (5) into Eq. (3) produces, for k e £3, 


T*+ 2T«+ 2 p*'Xf*'+c*X 2 i ki --J7 (.!*•**) -« ft Xc fc -h fc 


— f p X p dm = 0 (6) 

where either an overdot or the presuperscript i denotes an inertial reference frame 
for time differentiation of a vector* 

The integral in Eq. (6) cannot be evaluated explicitly without adopting a specific 
mathematical model of the flexible appendage; nor can one go beyond the integral 
representation of the substructure inertia dyadic J k and the vector c*, which vector 
establishes the shift of the substructure mass center c k relative to body-fixed point 
o k due to appendage deformation. Without commitment to a particular appendage 
model, one can accept (dropping subscripts within integrands) 

J = f (P • PU ~ pp) dm (7) 

and 



If tlie deformations u are to be assumed “small,” the term — X c k in Eq. (6) 
becomes negligible and can be ignored. 

Since it is possible to manipulate Eqs. (2) and (6) to eliminate redundant varia- 
bles and kinematic constraint forces and torques without making any restrictive 
assumptions concerning the appendage model or the size of its deformations, no 
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inhibiting assumptions will be imposed until specific cases are considered at the 
conclusion of this section. 

The immediate objective is the extraction from the 6n + 6 equations given by 
Eqs. (2) and (6) of 6 + n equations in the 6 + n unknowns established by the vec- 
tors X and «° and the scalars y,, • • • , y„ which define the relative rotations of con- 
tiguous rigid bodies. We must expect these equations also to involve unknown 
deformation variables for the appendages and rate variables for the rotors, but we 
must eliminate all interaction forces (typified by f ki ) and all interaction torques 
(typified by T'-'-*) except for those having the direction of the corresponding hinge 
axes, since these hinge-axis torques are known functions of hinge rotation or some 
more general control law. We must also eliminate the unknown position vectors 
typified by p r , replacing them by explicit functions of the hinge angles y u ■ • • ,y n 
and the system geometry, and eliminate all a> k for k e ff' in favor of terms involving 
Ci>° and y r for r € 

Simply by summing all equations defined by Eq. (2), we obtain the familiar result 

f = cynx ( 9 ) 

where Qii is the total system mass and F is the resultant of all external forces F r for 
Tcff}. Substituting from Eq. (9) back into Eq. (2) and rearranging produces 

F r + 2 f" - tyyir (F/(»? + p”) - o 

stfBv 


or, for r e£8, 


2 f" = -F r + C Vftr (F/Q 71 + P r ) (10) 


Summing over the branch set (see Def. 44) provides 

f’* =2 2 f rl = “2 [F r — Q?lr (F/042 + r)] 
so that, by Newton’s third law, 

f ki = -V* = 2 [F r — 0 Hr (F/Q41 + P r )] (11) 

r&kj 

Substituting Eqs. (10) and (11) into Eq. (6) eliminates all interaction forces from 
the rotational equations, providing 


2 [F r - Qvir (f /cm + p')] 


+ c fc x [-F k +cm k {F/m + p*)] - ■^(j‘*«») - qn& x <* 

pXpdm = 0 (kc£B) (12) 

J ft * 
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Equation (12) simplifies when written in terms of the vectors found in Defs. 33-36 
in the preceding section, due to the identity, for k 

2 p*' X 2 [F r - Wlr (F '/CM + pO] = 2 L*' X [F r - 0 9l r (F/QVl + p r )] 

fdflt r*&ki 

= 2 L fcr X F' + D kk X F - 2 L* r X (MrF 
= 2 (L* r + D**) XP-2L tr X QM r p r 

«£B 

= 2 (D kT X F r - L fcr X OWrPO (13) 


Substituting Eq. (13) into Eq. (12) simplifies its appearance somewhat, but there 
remains the problem of eliminating the unknown kinematic constraint torques 
which are present within the interaction torques typified by T*K By summing over 
all n + 1 vector equations defined by Eq. (12), we can, by virtue of Newton’s third 
law, obtain one vector equation involving no interaction torques at all; this summa- 
tion gives 


+ 2 (D* r X F r — L*' X 0?2 r p') + c* X [9H*(p* + F/OW) - F*] 
t «£B 

kt& 

- (J* * go*) - Oi?A-c f; X c k - h k - I pXp dm) = 0 (14) 


Equation (14), like Eq. (9), is free of kinematic constraint forces and torques, so 
that these two vector equations can be preserved in the final set that is our objec- 
tive. The remaining n scalar equations can be obtained by summing equations 
obtained from Eq. (12) over branch sets in order to isolate interaction torques such 
as T w , and then by dot-multiplying each expression by the unit vector parallel to 
the corresponding hinge axis. In particular, if we introduce the network elements 
found in Def. 32 of the previous section, we can sum over the branch set and 
get, for 3, 


T sy « + V'' ii* + 2 (r> fcr XT- L* r x Q?? r p r ) 

£—J 1 


Jtc#!*,* 


‘d 


+ c* x [cm t (p fc + f m) - 1*] - (j* • ta») - w xc‘-h‘ 


~ ftj » x r dm } ~ 0 

J fit, * 


(15) 


By Def. 32, N s < s, so that the labeled point on the hinge axis common to 6, and 
is, by Def. 11, called /«,; the unit vector parallel to the hinge axis is, by Def. 22, 
called g s ; and the interaction torque component of T 8 - v * parallel to g s is, by Def. 27, 
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Fig. 3. System geometry 


called t„. Hence the dot product of Eq. (15) with g E provides the scalar equation, 
for se£8. 


+ g 8 ' T + 2 (D*' XF'- L*' X Qn r p r ) + c* X [Qm (p* + F/0») - F*] 

^ V TtfPi 

— (J s * *>>*) — Q^fcC* Xc‘-h*--^j' pXp dm) = 0 (16) 




Equations (9), (14), and (16) provide the required 6 + n scalar equations free of 
kinematic constraint forces and torques, but they are not yet in final vector-dyadic 
form. The variable vectors typified by p r must first be expressed in terms of the 
system geometry, the deformation variables, and a subset of the 6 + n scalars 
required to define X, <o“, and the n angles of relative rotation at the n hinges. As a 
first step in this direction, we define & rj as the set of indices of bodies lying on the 
direct path between and A,-, and then (referring to Fig. 3) we can write 


p r - p 1 = c* + U r + 2 (L" - L*0 - I/' - c r (17) 

But L" = L rr = 0, by Def. 33, and for any index s in the set £8 — 6 r j the sum 
= 0, so that Eq. (18) can be written more simply as 

P r - P* = 2 (L* r - L") + c* - c r (18) 

st8$ 

Multiplying Eq. (18) by OWj/OVl and summing over all /e£8 produces 




ZZw< l, '- lw > + Z 




o % 

Oft 


(c> - c r ) 


je£3 


But 


y^Qw/p'=o 

m 
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by definition of CM, and 


so we have 



o % = cm 




sea 3 

Reversing the summation sequence in the first term, changing the index symbol in 
the final summation, and substituting the vectors introduced in Defs. 34 and 36, we 
find simply that, for re £8, 


-£( 


h sT + 




(19) 


Equations (14) and (16) contain not p r itself but combinations which, with Eq. (19), 
can be written as 


C* X X + |jp r ) “ *] 

rtSS 

and, in terms employing the symbols of Defs. 33-37, 

-SL*' X (ftlrp ' = — JjL*' X c tor r^(» ar + ^ fc ‘ S ) - ^ 


( 20 ) 


- -2 L* X CMr Td*' + s i> sr + - c r l 

- I v** J 

= 2 (D s * - D* r ) x [ <yn r i> kT + ( 2 Qn& - w&) ] 

re® L m «£8 J 

-2 2 % f L lf XD ,r - 2 012 ,L*'XD" 

tefB - k sefB-k-r rtfB—k 

..i/7 

= D** X 2 Q nXP* - -7 -(SD l, X 0M r D fcr ) +D»X2 0J2 8 c* 

f mB 


- D» X 2 On# - ( 2 D*9Mr) xVJ 

reffi re® it—/ W 


012,.. s 

c s 


+ 2D” cr X0llrC r - 2 2 QlJrL*” x D sr 

r«£B re® k «£8- fc— r 

+ 2 2 070^ x D rs 

r«£B-fc sefB-r 

= o - ^ [ 2 o>? r D* r x (t^ x d*') 1 + o-o + 2D»-x oj? r c r 


dt 


rt® 




“2 2 WxDn 2 2 

S€®-*-r rt®-k «€®-r 
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L d 

= - -T 7 [{ 2 (D kr • D fcr U - D^D^)} • *>*] 

(It fg£3 

+ 2 QWrt>* r x c r +2 2 QWs (l^ r - L &8 ) X b r * 

r*£*> r*&—k sffiS—r 

or 

l d 

-2 L fcr x Oft r p' — j- [K* • ca*] + 2 (MrD kr X ‘d- 

rt£8 

+ 2 2 QW. (L* r ~ L fc ») X D r « (21) 

Kffl -k «/® 

In the development leading to Eq. (21), it was recognized that 

2 Qw£> r8 = 0 

s*m 

(reflecting the significance of the barycenter as the mass center of the undeformed 
augmented body, which consists of the substructure augmented at each connection 
point by a particle having the mass of the corresponding branch set). It was also 
noted that 

-2 2 Q^rL* r x D*' - -2 2 QMrL*' x i> 

k— r se$~k re£ft—k~s 

= -2 2 %L**XD r ' 

re®— k r 

where the second step involves simply relabeling indices. 

Finally, we can recognize in Eq. (21) that the quantity U ir - L* s is zero for any 
index s corresponding to a body which lies anywhere on the branch which begins 
with ^ V);r and includes i T , and for any other index the quantity D" is also D r *. Thus, 
in Eq. (21), D™ can be replaced by D rfc , with the consequence 

‘d 

~ 2 L kT x o 7lrP T = ~ -r (K k • u*) + 2 QM, D kr X c' 

+ 2 2 Q"h s (L kr - W) x fr* 

«©- k 

= - -j7 (K ,f • «*) + 2 QW,D fc, ‘ X c v 

at «ffl 

+ 2 [((MS' - 2 WhL ks ) X D Wt ] 

r*£8~k 

or, with Defs, 34 and 36, 

l d 

- 2 L fcr X (yn T p r = - -T7 (K* • <■>*) + 2 9ftrD* r X c v + 2 Q4?D fcr X D 1 * (22) 

By substituting Eqs. (20) and (22) into Eqs. (14) and (16), we can obtain the 
desired vector-dyadic equations. In this substitution, it becomes apparent that K* 
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and J fc always appear in combination, suggesting the introduction of the new dyadic 
(see Def, 38). Then the new versions of Eqs. (14) and (16) become 


T; jT* + SD*' XF' + C^X (^ F - F k ) + D fcr X c r 

+ QW*c k Xs(dH Ttr'c') - -i (•*■ W*} + 2 Q«D* r X D A - h* 

\ \Yn J at 

— j p X p dm | = 0 (23) 

and, with judicious use of the path elements found in Def. 29, 


T, + gs * 2 (THSD* r XF' + c fc X(^F 

kt& l re® \ ' ll 


/Wh 


- F*) 


+ 2 ^rD* f xc r + cm k <* x ^ (ir* + ^ fcV ) - ^ (•“ * w *) 

+ s W'XD*-ii'-^) p X p dm | — 0 (24) 

«®-* at ) 

Equations (23) and (24) can be restated in more useful form by expanding terms 
involving time derivatives relative to inertial space to obtain time derivatives rela- 
tive to reference frames established by individual substructures, plus additional 
undifferentiated terms. In particular, noting that D 7 * is fixed in 4 r , we can substitute 


D r * = &x T> rk + <0 r X (co r X D rk ) 


(25) 


so that in terms of Def. 39, we have 

2 Q»tf f xfr fc =-2 71 2 D fcr X [co r X (co r X D r *)] (26) 

n&-k 

The parallel expression for c r becomes 

c r = °c r 4- 2co r X c r + X c r + o> r X (co r X c r ) (27) 

with open circles representing time differentiation in the reference frame estab- 
lished by br- 
other time differentiations in inertial space expand as 


o 



h fc = hH X h fc 

(28) 


. o , 


£<**"*> 

= • <0 & -j- <0* ' X • CD & 

(29) 
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and 


Vi 

~di 


/ pXp dm = / 


p X p am 4- co* X 



dm 


(30) 


with the open circle indicating time differentiation in the reference frame estab- 
lished by the local substructure (here 6k). 


By combining Eqs. (23-30), we can obtain the vector-dyadic equations of rota- 
tion in the form 


and 


where 


2W* = 0 


r, + g 4 * 2^w* = o 

kt& 


(31) 

(32) 


W* * T* + 2D* r X F r + «*x(^ F -F*) 


+ 2 0 VlXP r X + 2w r X c r + (o' X c r + co r X (a> r X c r )] 

rt0R 

+ tyftkC* X2[d'X D rk + CO' X («' X D r *)] 

+ Ofec* X ^ ^ ^ [c >r + 2to r Xc' + u r Xc' + « r X (w r X c r )] — <!>** • m* 

rtm 

- 2 (**' * C 0 r ) + <Yn 2 D fcr X [w r X (co' X D r *)] - co* X <&** • CO* 

re®— Jk 

- h* - to* X h* - O** -to* - / p X p’clm - co* X f (pXp)dm (33) 

“'"k 

To obtain the vector-dyadic equations in their final form, involving the 6 n 
unknowns in X, w°, and yi, ‘ ■ • , y n as well as the rotor momentum variables and 
the deformation variables, we must substitute in Eq. (33) the expressions 


and 


<0* = 40° + 2 Crkir g' 

nff 


(34) 


co* = <b° + 2 [y«-g r + x yrg r + 2 e,ri»g s x gl (35) 

re# «# 

where e r * are the path elements found in Def. 29, and Defs* 9, 22, and 23 are also 
employed. 


In combination, Eqs, (9) and (31-35) provide the final form of the vector— dyadic 
equations of vehicle translation and substructure rotation. In addition to these 
equations, for each rotor in the system an additional scalar unknown is introduced, 
and one more differential equation is required* If a rigid axisymmetric rotor with 
spin axis inertia spins relative to body at the angular rate fas about an axis 
parallel to the unit vector h k * fixed in 4 ky and if r % 9 is the magnitude of the b k8 com- 
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ponent of the torque applied to the rotor, then the required additional equation of 
motion is 

(36) 

If there are N H rotors in <4, and N B rotors in the entire system, then Eq. (36) applies 
for s = 1, ■ • ■ ,N k , with k e £6, so that Eq. (36) contributes N H scalar equations. 
Moreover, the variables in equations typified by Eq. (36) are related to those in 
Eq. (33) by 

h* = 2 J%ks$kab ka (37) 

8 ~ 1 

The scalars typified by r£* must of course be given, either explicitly or implicitly 
in the form of additional differential equations representing control laws. The same 
is true of r 3j T* and F r in Eqs, (33) and (32). Once these issues are settled, however, 
there remains the problem of defining the nonrigid appendages mathematically 
and establishing their vibrational motion equations. This issue has been postponed 
until the end of this section, in order that we might consider several possibilities. 
The relative advantages of alternative appendage models are considered in detail 
in Ref. 6 and are only briefly summarized here. 

If the appendage a k is of extremely simple configuration (such as a uniform elastic 
beam), and the nominal angular velocity u h of the base to which it is attached is of 
a special class (such as zero, or parallel or orthogonal to the hypothesized beam), 
then it might become attractive to model the appendage as an elastic continuum, 
writing partial differential equations of vibration, seeking normal modes of vibra- 
tion, and truncating to a modest number of coordinates and a corresponding num- 
ber of ordinary scalar differential equations. One must then return to those terms 
in Eqs. (31-35) which depend on deformation (noting Eqs. 7 and 8) and evaluate 
those terms as functions of the modal deformation coordinates, ignoring second 
degree terms in deformation variables. 

In most cases, it will prove more feasible to adopt for an elastic appendage a 
finite-element model such as that described in Section I. The equations of small 
vibratory deformation can then be recorded from Ref. I 4 (see Eq. 64), and the 
appropriate coordinate transformation can be adopted from Ref. 1 (see Eq. 74). 
The vector c* in Eq. (8) is then replaced by c as recorded in Ref. 1, Eq. (43), with 
e = 0 and summation extending only over the fcth substructure. Although the 
dyadic J /l in Eq. (7) and the integral 

if 

J «k 

are not provided in Ref. 1, which treats the distributed-mass finite-element model 
with nodal bodies, they are available in Ref. 5 for the special case with all ma^s 
concentrated in the nodal bodies; see Ref. 5, Eq. (126) for J* and Eq. (127) for 
and see Eq, (112) or Eq. (114) for 


fd 

dt 


L 


pXp dm 


4 Note that the vector called X in Ref. 1 must be interpreted for the rth appendage of the present 
paper as X + p r , with p ? expanded as in Eq. (19). 
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recognizing in each case that the superscript k indicates restriction to the kth sub- 
structure, and setting Sl a and | as found in Ref. 5 identically to zero. In each case 
the indicated expressions involve deformation variables which must be subjected 
to coordinate transformations and truncations, as indicated in Refs. 1 and 5. 

If it pleases the analyst, he can conceal the complexities of modeling the elastic 
appendage by simply adopting a symbolic model in terms of mode shapes and 
natural frequencies. Then he can rather easily record the equations of vibration 
(probably using Lagrange's equations with modal coordinates for generalized 
coordinates), and he can evaluate the deformation-dependent terms in Eqs. (31-33) 
in terms of unspecified mode shapes and frequencies. Of course he must face the 
problems of adopting an appendage model and determining its modal character- 
istics when he wants to face any real problem. 

VI. Matrix Equations 

This section illustrates the transition from the preceding vector-dyadic equa- 
tions to an alternative matrix form which is more explicit and therefore more suita- 
ble for programming for digital computer numerical integration. To obtain explicit 
results, we must make a commitment to a particular appendage model, sacrificing 
some of the generality of the vector-dyadic equations. In what follows, we confine 
attention to the special case of the finite element appendage model in Ref. 1 for 
which all mass is concentrated in the n nodal bodies of the system (with no dis- 
tributed mass for the internodal elastic elements). All deformations from a nomi- 
nal appendage state are assumed arbitrarily small, so that terms above the first 
degree in these deformations and the corresponding deformation rates can be 
neglected. Moreover, we assume that only the reference body contains rotors, 
and these number three, establishing an orthogonal triad whose axes parallel b°, b°, 
and b? v The restrictions adopted in this section reduce the number of terms appear- 
ing in the final equations of motion somewhat, but (as shown in Ref. 1) they do not 
change the structure of the equations, so that the special case treated here remains 
representative of the general case of small-vibration appendages. (Large-amplitude 
appendage vibrations, which are accommodated in the preceding vector-dyadic 
equations, would produce matrix equations quite different from those to be 
recorded here.) 

The vector-dyadic equations to be written in matrix form are Eqs. (9) and (31-35), 
in combination with Eq. (64) of Ref. 1, with various results substituted from Refs. 1 
and 5. 

Equation (9) may be rewritten in terms of Defs. 8, 15, and 28 as 

{i} T C T F = {i} T (yftX 


which implies the matrix equation 


F = mCX (38a) 

The time history of the direction cosine matrix C is established by the Euler— Poisson 
kinematic equation (see Ref. 5, Eq. (44), for tilde convention), 

C = (38b) 
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Skipping next to the rotor equations and noting the restriction to the three orthog- 
onal rotors in (so that k = 0 only, and b 05 = b° for s = 1, 2 ? 3), we have from 
Eq. (36) the tliree scalar equations 

Tm =Jfos('£os + «J) 


or the implied matrix equation 




(39) 


where ^ = [^ 0 i ^02 ^os] and ^is a diagonal matrix with diagonal elements Jf 01 , $> 2 , 
and 

The vector rotational equations are easily written in matrix form with the 
definition 

W* = W* ■ {b*} 

since then, in terms of the direction cosine matrices in Def. 31, we can write 
Eqs. (31) and (32) as 


or simply 
and 


or 


2 {h h ) T W k = {b°} r 2 C ok W k = 0 
urn mb 

2 C°*W* = 0 

MB 


r, + g aT {b'}^2^{b fc } r W k =0 

M? 


r. + g- T 2^c-w = 0 


(40) 


(41) 


Equations (40) and (41) require W\ which is available from Eq. (33) when writ- 
ten in the form 


{ b h } T W k = {b k ) T T k + 2 {b k } T D kr X {b r } T F T 

T£& 

+ (^{bTF- (bT'F 1 ) 

+ 2 Qftr {b*} r D kr X {b r } 7 ' [c r + 2z r c r + Z r c r + t> r M r c r ] 

-I- Ofe { b k } T c k X 2 {b r } T [£ r D rfc + - (V} T ^V 

re® 

+ 2 (OK [D r * r {b r } • {b*} r D kT {b k } r U {b*} 

n£S-k 

- {h r y D rk D krT {b fc >] • {bT" r } 

+ cm 2 [ (b fc } T D* r X {bfSVD'*] - (b fc } r V 

niS-k 

- {h*} 2, h k - {b*} r 3W ! - {b*} r f ppdm - {b fc } r <3* / pp cim 

1/ •/ 
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since then 


W*= (b^-fbTW* 

= T k + 2 D kr C kT F r + C' CV F - F k \ 

*a \on / 

+ 2 Q‘ftrD kr C kr (c r + 2s T c r + Z r c r + s r s r c r ) 

reffi 

+ tylgp 2 C kr (Z r D ,k + «SVD rJt ) ~ $**&* - &*J> 

+ Q71 2 [(D^C^D^t/ - C kT D rk D krT ) C kr i> r 

rdB-k 

+ D fcr C 4 ’o> r sr r D rft ] 

— asMV — /i* — m k h k — / pp dm — a' 1 I pp dm (42) 

Even with the substitution of Eq. (42) into Eqs. (40) and (41), the matrix equa- 
tions are not in their final form, since the variables a> k , m k , h k , c k , #** the direction 
cosine matrices such as C kr , and the integrals over a k are not yet sufficiently explicit. 
These deficiencies are rectified by substituting matrix versions of Eqs. (34), (35), 
(37), (8), and (7), together with new expressions for C' !T and expressions for the 
integrals obtained from Ref. 5. Specifically, from Eqs. (34) and (35), we find 


= {b*} *w* = (b fc ) • [<bo}*V + 2 ^{bTg'l 

f«£P 


or 

= C fc V + 2 erky r C kr g' (43) 

r tff 

and 


= {b*} • c 2 » fc = {b fc } • (b“} ! ’ to® 


or 


+ 2 [y, {b*} • {b'} T g r + {b*} ♦ {h 0 } 7 <»° X {fry g' 

rtd> 

+ {b‘}‘S^{b , } r g , X{bTg1 

Sf£P 


& k = C*°<o° + 2 «* [yrC kr g r + C*VC°'g' + 2 «.rV *C k ‘g‘C° r g r ] (44) 

rtff $tff 

For the special case at hand, Eq. (37) yields h k == 0 for k ^ 0 and 

b° = {b°} • h° =Jf<p (45) 

as these symbols are defined following Eq. (39), 
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For the nodal body system to which this section is limited, Eq. (8) becomes a 
special case of Eq. (45) of Ref. 1 or a variant of Eq. (60) of Ref. 5. Specifically, 
we have 



8 — 1 


where appendage has been idealized as n k nodal bodies interconnected by mass- 

A . 

less elastic structure, with m 8 the mass of nodal body s and u* = [t#{ uj u 8 3 ] T given 

by u s 4 (b fc ) *u s with u 8 the displacement of the body s relative to from the 
position occupied in the nominal state. 

Recalling from Defs. (38), (37), and (24) that 

- K* 4- J* (47) 

and noting that K k is constant and J k is available from Eq, (7) as the variable matrix 

J h — I (p T pU — pp T ) dm 
•h j; 

one can, for the presently adopted appendage model, evaluate } k as (see Ref. 5, 
Eq. (126)) 


p = P-+ 2 {m. [2 (R k + r s ) T u*U - (R* + r°) u* T - W (R k + r*) r ] + fi‘l‘ 


I*p*} 

(48) 


where ] k is the nominal (constant) value of J k , and I s is the constant inertia matrix 
of the sth nodal body in its own body-fixed vector basis {n s } A , where, in the nominal 
state, {n 8 }* = (b'} . Here we use u s for (u 8 )\ as in Def. 21. 


In combination, Eqs. (47) and (48) provide 


<*>»* = 2 {m, [2 (R ?; + i s ) T u‘ U-(R k + r*) u sT - u‘ ( R k + r*) T ] + p*I* - I s R s } (49) 

8 = 1 


An expression for the direction cosine matrix C k2i k (where k e £75 and N k is a net- 
work element from Def. 32) may be found from Ref. 7 to be 

C kN k = U cos y* — g* sin yii + g*g* r (1 — cos y k ) (50) 

and, by transposition, 

CV = U cos y k + g* sin y* + g fc g* r (1 - cos y k ) (51) 

For arbitral r, / €£3, the direction cosine matrix C r > can be constructed as a product 
of direction cosine matrices relating contiguous bodies (as given by Eqs. 50 and 51). 
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The required product can be written symbolically as 

Qri = II IVp™* (52a) 

where the symbol IJJV is understood to imply the following algorithm: 

A 

(1) Define p = N rj and construct C rp from Eq, (50) if r > p and from Eq. (51) 
if r < p. 

A 

(2) Define q = N P j and construct C™ from Eq, (50) if p > q and from Eq. (51) 
if p < q. 

A 

(3) Proceed until an integer u emerges such that / — N Up finally constructing C n j 
from Eq, (50) if u > / and from Eq. (51) if u < 

(4) Multiply the matrices obtained in the sequence 


C ri = CC* * * * C ui (52b) 

Finally, Eq, (42) requires more explicit expressions for the integrals over the 
appendage For the special case of elastic appendages idealized as intercon- 
nected nodal bodies, the appropriate expressions can be found in vector form in 
Eq. (114) of Ref 5. The corresponding matrix equations are 

R* (2 + w* 2 — 2 (ZPr 8 )" m 8 ii* 

2 r* (' m a u * + — 2 I s 0 a + « fc /? fl ) 


— jf pp dm — I pp dm — 


- 2 1'P* + 2 


~ ”2 (R fc + 7 s ) - 7» k 2 + ? a ) m s ti 8 


“ 2 U s (/?* + !•/§* + (53) 


where 2 implies summation over the range from s ~ 1 to s = n k , and the super- 
script k has been dropped from nodal body variables in the fcth appendage (such 
as u\ which replaces (u 8 )*). In developing Eq, (53), the identity (3fc)~ =ab~ba 
has been utilized. 


In summary, in order to obtain the explicit matrix equations defining substruc- 
ture rotations, one must combine Eqs. (42-53) into Eqs. (40) and (41). The result is 
a set of n + 3 scalar equations, which can be collected together with the three rotor 
equations (Eq. 39) and the three vehicle translation equations (Eq. 38a) and their 
kinematic equations (Eq. 38b) to comprise a set which is complete except for con- 
trol system equations and the equations of vibration of each of the n + 1 append- 
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ages. The latter set can be obtained for a distributed-mass finite-element model 
from Eq. (64) of Ref. 1, and for the special case of the concentrated-mass (nodal 
body) finite-element model either from Eq. (64) of Ref. 1 or from Eq. (95) of Ref. 5 
(correcting the last algebraic sign within the braces on the right side of Eq. (95) by 
changing — to +), recognizing that what is recorded as ©X in those equations 
must be replaced by the inertial acceleration of the mass center of the correspond- 
ing substructure in the local vector basis, which for substructure is given by 

CX + {b*} • P* = CX+ {b*} T XX r ) - **] 

= CX — (c* + + 2S*c* + «Vc fc ) 

+ 2 C*- [(S’- + v#)(d* + H cr ) + yfif ( 6r + 2sr ^ r )] 

(54) 

where Eq. (19) has been utilized for the expansion of p fe . 


Typically the number of nodes n k in a single finite-element model of an elastic 
appendage is so great that the proposed combination of nodal body vibration equa- 
tions with substructure rotation equations, vehicle translation equations, rotor 
equations, and control system equations would be prohibitive in size for numerical 
integrations. Thus it is to be understood at the outset that the nodal body vibration 
equations will provide the basis for a transformation to distributed or modal coordi- 
nates for appendage deformations, and that most of the modal coordinates will be 
deleted from consideration by truncating the matrix of deformation variables. By 
means of this transformation, all variables labeled u* and in Eqs, (46-53) are 
replaced by appropriate combinations of new modal deformation variables. Thus 
the equations actually to be programmed for digital computer numerical integra- 
tion are not precisely those obtained from Eqs. (38—53) of this report plus Eq. (64) 
of Ref. 1 (or Eq. 95 of Ref. 5), but instead the transformed and truncated versions 
of these equations. Because questions of transformation, truncation, and synthesis 
of the resulting composite matrix equations have been treated extensively in Refs. 1, 
5, and 8, these matters will not be considered further here. 

VII. Summary 

By combining Eqs. (9), (31), (32), and (36) with Eq. (64) of Ref. 1, one can obtain 
a generic statement of a minimum dimension set of equations of motion of a system 
of n + 1 rigid bodies interconnected by n line hinges, with a set of axisymmetric 
rotors and a distributed-mass finite-element model of an elastic appendage attached 
to each rigid body. In order for these equations to stand alone as a complete formu- 
lation of the problem, one must substitute the auxiliary equations (33-35) and (37), 
as well as certain specified equations from Ref. 5, into Eqs. (31) and (32), and 
Eq. (64) of Ref. 1 must be composed from the underlying equations of that ref- 
erence (Eqs. 46, 53, and 62). Moreover, the total system of equations must for 
practical utilization be written in matrix form and the appendage deformations 
subjected to modal coordinate transformations and truncations, as these procedures 
are described in Refs. 1 and 5 and illustrated in part for a special case in the preced- 
ing section of this report. 
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Although one must not underestimate the labors of proceeding from the vector- 
dyadic equations in this report to a generic computer program for integrating these 
equations, this does appear to be a feasible task. Once the program is completed, 
it will have sufficient generality to encompass several spacecraft simulations for 
which specific numerical integration computer programs have been written in the 
past (Refs. 8-10). Experience with these programs provides reasonable assurance 
that the generic program proposed here will find application as a practical tool for 
the simulation of complex modern spacecraft. 
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